Exact results for a Kondo impurity with electronic correlations 
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An one dimensional model of a Kondo impurity with arbitrary spin in the presence of electronic 
correlations is solved exactly. The Bethe ansatz equations are obtained for the case of quadratic 
dispersion. The ground state is studied in the thermodynamic limit, and the magnetic susceptibility 
at zero temperature is calculated. 
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It is known that the study of exact solutions is use- 
ful to the understanding of non-perturbative effects in 
the strongly correlated electronic systems. The exact 
solution of the one dimensional Kondo model with lin- 
earized dispersion and no correlations between electrons 
was found in |^,^. In this paper, we study the problem 
of a Kondo impurity with arbitrary spin in the presence 
of the electronic correlations. We show that the model is 
exactly solvable at some value of electron-impurity cou- 
pling. Using Bethe ansatz, we solve the secular equations 
for the spectrum with periodic boundary condition. The 
ground state as well as the magnetic susceptibility at zero 
temperature are calculated. 

We focus on a system described by the Hamiltonian 

k 

+ X! Wf^lkl +k2,k3 -|-k4)Ck^^Ck3aCk2aCkia, 
ki,k2,k3,k4 

where Cka annihilates an electron with momentum k and 
spin component a. The electrons are coupled by both 
spin and charge interactions to a localized impurity, 

Hi = ./*:(0)S,fc*fc(0) • 8° + y*:(0)*„(0), 

where the field 'i'a is the Fourier transform of Cka, S° 
is the spin of the impurity and S is the spin of the elec- 
trons in the band. If we restrict our attention to the 
low-energy regime, namely, k close to the Fermi surface, 
we obtain the (linear dispersion) one-dimensional Kondo 
Hamiltonian with a J-function interaction between elec- 
trons. The coupling between the electrons docs not con- 
tribute to the scattering matrix (S-matrix) in the case of 
linear dispersion. However, once correlations are intro- 
duced, the scattering matrix is determined uniquely from 
the Schrodinger equation and the anti-symmetry of the 
fermionic wave function. 

We consider the case of quadratic dispersion, i.e., 
e(fc) = (in units of h and of the electron mass). 



We show that the present model can be solved exactly 
by means of Bethe ansatz with periodic boundary con- 
ditions. In order to have a periodic Hamiltonian in the 
Hilbert space of TV electrons, we define the following an- 
ticommutation relations for the fields, 

{*(a:), '^'^y)} = 5{x - y - nL). 



Then the Hamiltonian in the Hilbert space takes the form 
H 
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^ (JS^' • S° + V)6{x.j - nL). (1) 
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The Hamiltonian operator (|l|) is invariant under arbi- 
trary permutation. Although it does not have transla- 
tional symmetry due to the presence of the impurity, (|l|) 
is invariant under a translation of finite distance L. Let 
{Pq} be the set of hyperplanes Pq := {x|(a|x) — 0}, 
where the notations x := {xi,X2, • • • , xn), a = a — Cj, Ci, 
and {a\x) := '^f^i{ct)iXi are used. In the domain \ 
{Pq}, eq.([^) becomes an iV-dimensional Laplace opera- 
tor, so iV-dimensional plane waves solve the Schrodinger 
equation. As those hyperplanes partition the domain into 
finitely many regions, the wave function must be a piece- 
wise continuous function. These regions can be labeled 
by both an element of the permutation group Sn and 
a number indicating the ordering of the impurity and 
the electrons, namely, C{Q^'^^) := {x\xqi < ■ ■ ■ < xqi < 
< • • • < xqat}, here i = 0, 1, 2, • • • , and Q^^^ means 
< xqi < ■ ■ ■ < xqn. In these notations we can write 
the wave function in the Bethe ansatz form. 
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where a := (agi, aQ2, • ■ • , o-qn), o-j denotes the spin com- 
ponent of the jth particle; Pk represents the image of 
a given k := (/ci, A:2, • • • , /cat) by a mapping P e Sn] 

{Pkm^Ef=iiPkUQ^h- 

The wave function for the electrons must be anti- 
symmetric under any permutation of both the coordi- 
nates and the spin components of the electrons. Since 
any element of the permutation group Sn can be ex- 
pressed as a product of the neighboring interchanges, : 
(• • • , Zj, Zj^i, • • •) 1-^ (■ ■ ■ , Zj+i, Zj, • • •), where zj denotes 
either Xj,aj or their image by a mapping of Sn- The 
requirement of anti-symmetry is {a■^^p)a{x) = —ipaix). 
Using the identity {Pk\a^x) = {a^ Pk\x), we obtain 



(3) 



where the spin labels are omitted and pQi'Q(i+^^ is the 
spinor representation of the permutation . Explicitly, 
we have V'^^'^^^+^'> (1 + S'^^' • S°)/2. The function 
terms in the Hamiltonian (|l|) contribute boundary con- 
ditions at hyperplanes {Pa} |3[^, i.e., the discontinuity 
of the derivatives of the wave function along the normal 
of the hyperplanes: 

^lirn^ {a ■ V [ipa{x(a) + ea) - il^aix^a) - ea)] 

-U [lpa{x^a) + ea) + iJa{x(a) 



where V := J2f=i'^ji'^/'^^j)' ^ 



(a) 



- } = 0, (4) 

€ Pq, L/ = M for 
a = e^- Cj and U ^ JS^ -8° + V ior a ^ e^. 

As the regions C(Q(*)) and C{a^Q^'^) {j ^ i) are ad- 
jacent to each other, and so are C(Q'*') and C{Q^^~^'>), 
the respective boundary conditions (Q) relate the wave 
functions on neighboring regions. After writing out the 
boundary conditions (0) and using (^, we find that: 
A(P;CTJg«) = (Pfc)j+i]A(P; Q^) for 

j ^ i, and A(P;g(*^i)) = S'^^'0[(Pfc),]A(P; g«). The 
S-matrix S'^-''^ of the electron-impurity depends on the 
coupling constants J and V, and the S-matrix of electrons 
takes the same form as that of the (5-function fcrmion gas 
ii- 

Both the electron-electron and electron-impurity S- 
matrices relate the coefhcients A between distinct regions 
in the configuration space of N electrons. Using (||) and 
the S-matrices, we find that the coefficients A in any 
region are determined up to an overall factor by the S- 
matrices, ^-^^-QO+i) := _7?QxQW+i)5Q^,Q0+i)^ 
i.e., A(a*P;Q«) = 5Q(*+i).o ^'S'.Ql^+i) 
A(P;gW), and A(cr^P;g«) = 5'3J'Q0'+i)yl(P; gW) for 
j i- 

The Yang-Baxter equations arising from both 
A(CTJCTJ+VJP;gW) = A(CTJ+VV^'+iP;gW) and 
A(P;CTJcr^+VjgW) = A(P;crJ+VJcrJ'+ig«) for i ^ 
j,i -f 1 are fulfilled identically. However, the others 
for i = j,j + 1, namely, 5'Q^Q(»+i)5Q^05Q(»+i),o = 
5Q(»-i-i),05Q»,05.Q»,Q(»-i-i) g^j.g satisfied only when J = 

and V arbitrary, or J = —2u and V = —{s + l)u or su 



for an impurity of spin-s. The former corresponds to a 
non-magnetic impurity Q . In the following we will focus 
on the later case. The electron-impurity S-matrix reads 



-te[{Pk),] 



i{Pk\ 



(5) 



where T^'^^ = (1 + S*?^ • S°)/2, 6{x) = 2tan-i(a;/u(s -|- 
1/2)) + TT ioYV = -{s + l)u and e{x) = for F = su. 
These two cases correspond respectively to the anti- 
parallel and parallel spin-coupling between electron and 
impurity. 

Now we proceed to determine the secular equations 
for the spectrum by considering periodic boundary con- 
dition. If a; is a point in the region C(g*^*''), then 
x' = (xi, • • • , xqi -I- L, • • • , xn) is a point in the region 
C(7g(*-i)) with 7 = cr^-V^-2...^2^i^ ^Yic peri- 
odic boundary condition imposes a relation between the 
wave functions defined on C(g(')) and C(7g('-i)). Writ- 
ing out this relation in terms of (H), we find that the pe- 
riodic boundary conditions are guaranteed provided that 
yl(P;7g(*-i))e*(^''')i^ = A(P;g(^)). After applying the 
S-matrices successively, we obtain an eigenvalue equation 
in the spinor space: 



gQl,QN ^ ^ ^ gQl,Q(i+l) gQl,0 gQl,Qi . . . S'^^'^'^A{P- g^'^) 

^A(P;g(')). (6) 



Applying the quantum inverse scattering method (see 
[Q), we find the Bethe Ansatz equations: 
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Ay — + iu Ap — ius 4\ \^ — ki — m/2 
Ay - A^ - iu A^ -f ius ^J-^ - ki + iu/2 



(7) 



Taking the logarithm of (^, we obtain a set of coupled 
transcendental equations. 
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5]ei(Ay-/c,) + 7?e,+i(fc,) 

AT 

^ei(A^- Ay) - -27rJ^ + ^ei(A^-A:0 + e,,(A^), (8) 
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where 6/3(2;) := 2tan~^(x//3M), — vr < tan~^(a;) < tt; 
?7 = 1 for F = — (s + l)it and 77 = for V = su. The 
spin quantum number takes integer values or half- 
integer values according to whether A'^ — M is even or 
odd; The charge quantum number Ij takes integer values 
or half- integer values according to whether ill -I- 1 (or M) 
is even or odd for V = — (s -I- l)u (or for V = su). The 
last terms on both right hand sides of (P) represent the 
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contributions of the impurity. The phase shift due to the 
impurity affects the quantum numbers Ij and J^. 

We consider the thermodynamic Hmit. The ground 
state of the present model is a Fermi sea described by 
Poik) and <7o(A), where po{k) is the distribution function 
of charge with momentum k and (Jo{X) that of down spins 
with respect to the rapidity A. The subscript zero is used 
for the case of zero magnetic field. The distributions of 
the roots satisfy the following coupled integral equations, 



Poik) 



1 

2^ 

1 
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Ks+iik) + 



dX'(7o{X')Ki{k- X'), 



-B 



dX'ao{X')Ki{X~ X') 



D 



dk'po{k')Ki{X- k') 
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where K/jIx) := I3u/{(3'^u'^ + a;^)7r. 
determined from the conditions 



(9) 



The B and D are 



J dX ao(X ) = — , 
J dk'po{k') = — . 



The energy of the ground state can be calculated once 
po{k) is known. In the absence of magnetic field B = oo, 
we can write (j^) in a closed form by Fourier transform, 



27rpo(fc) 



47r 

u ' 



Do 



dk'po{k')Ri 



-Do 



2k - 2k' 



where 
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In the strong coupling limit u ^ 1, we obtain 
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where = (-l)'*[l/4 + Ei=i(-1)7(2^ - 1)7^1 for im- 
purity with s = integer and Ws = (— 1)'*^ 2 [In 2/27r + 
l)'/2/7r] for impurity with s = half-integer. 
From the explicit form ( |ll|) we can completely determine 
D, and the energy density of the ground state is 
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In the presence of a magnetic field, we add to the 
Hamiltonian (|l|) a Zeemann term, —hJ2j'=o ^i- 
add a term arising from the nonvanishing chemical po- 
tential. For h 0, B depends on the magnetic field. It is 
convenient to introduce the dressed energies of charge 
and spin, ep(fc), edX), which are defined by 



E 
T 



dkp^°\k)ep(k) 



rfAf7(")(A)e,(A)-/iy, (13) 



where p(0) = l/2Ti-r]K^^i{k)/L andcr('')(A) = Ks(k)/L. 
The equations (^) of Fredholm-type |^ for p{k) and a{X) 
are equivalent to the following integral equations for the 
dressed energies. 



e^(fc) =e(°)(fc)+ / dAe,(A')i^i(A'-fc) 



e,(A) =4°)(A)+ / dk'ep(k)Ki{X' ^k) 

J-D ^ 

dA£,(A')i^i(A'- A), (14) 

-B 

where ep"^ are the bare energies, namely ep°'' = k'^ /2+p — 
h/2, elj ^ — h. The solutions of ( [l^ ) define the pseudo- 
energy bands. The filling of all states with ep(fc) < and 
£ct(A) < corresponds the ground state. The Fermi levels 
of both fc-sea and A-sea are determined by the conditions 
£p{D) = and £a{B) = 0. 

Once the distribution ct(A) is solved, the magnetization 
per unit length follows as. 



M = 



' dXaBiX) + h^-h. 



(15) 



Clearly, A4 is not exactly zero when B = 00 and s 1/2 
due to the magnetization of the impurity. We will focus 
on the case of the weak magnetic field. Then the o'(A) for 
large A mainly contributes to the magnetization. With 
the assumption X ^ D, the second equation of be- 
comes 

ctb{X)+ f dX'IC{X-X')aB{X')^a^{X). (16) 

J\\'\>B 

where a new kernel /C and a distribution function (Too (A) 
have been introduced. They satisfy the following equa- 
tions 



/C(A - A') = Ki{X -X')+ / dA"i^i(A - A")/C(A" - A') 

■/ — 00 

<Joo{X) ^ a^"\X) + dX'KiiX- X')aoo{X'). (17) 



Shifting the origin, by introducing a'(A) = asiX + B) and 
neglecting the higher-order terms, we obtain the follow- 
ing Wiener-Hopf integral equation for ct(A). 
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cr{X)+ d\'IC{\-X')a{\') = croc{X + B). (18) 

Equation dl^ ) can be solved by the one-side Fourier 
transform |9|, and reads 

{V2^IC{uj) + l)a+ {uj) + a- (u) ^ aoo{Lo)e-"^^ (19) 

with 



(20) 



(V27r/C(w) + 1) 



2cosh(fw)' 



where ij^{lo) and are the analytical parts of the 

Fourier transform u{lo) on the upper and lower w-plane, 
respectively. The kernel ( pO| ) is regular in the strip 
— - < Iniw < -, and is factorized into /C+(w) and JC^(iS) 
such that 



(21) 



equation for £cr(A) from the second equation of (14). Us- 
ing the similar procedure, we solved the dependence of 
B on h. An explicit expression for the susceptibility at 
zero temperature ( p2| ) follows with dB/dh = —l/irh. 

In the above we obtained exact results for the Kondo 
impurity by taking into account the electronic correla- 
tions in the case of quadratic dispersion relation. We 
showed that the model is integrable when the correla- 
tion strength u is proportional to the J, the strengh of 
electron-impurity coupling. If we consider the approx- 
imation kj ~ ki for any j, I, the S matrix of electron- 
electron will be independent of u. Then the Yang-Baxter 
equation will give no relation between u and J. This 
makes it easy to understand the usual Kondo problem 
where the linear dispersion relation is adopted. 
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{V2nIC{uj) + 1) = IC+{lu) JC-{uj) 



where /C+ = {-izf^ e-'^T{l/2-iz) / with z = ww/27r 
and [AC^(a;)]^-^ — Ai+(— w). AC+(ijj) and K~{lo) are regu- 
lar and free of zeros respectively in the upper and lower 
half-planes. After seperating Yi~ (Lo)Ua:,(u))e~'^'^^ into a [1 
sum of two parts q+iiS) + q-(u)\ which are analytical in 
upper and lower half-planes respectively, we obtain the 
solution of ( p^ i5-+(w) = qj^{ix:) j K,'^ [lJ) . The magnetiza- 
tion at zero temperature is 



M 



^-(o) + ^(.-i). 



We infer the magnetic susceptibility as 



\ \ dB 



n=0 



■fs{n)e 



-(2n+l)f S 



(22) 



where /s(n) = exp{— i(n-|-i)[l — ln(m-t-ii)-f (2s— l)7r]}. 
The relation B = B(h) follows from ea{B) = 0. In the 
limit of strong coupling, we also obtain a Wiener-Hopf 
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